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Abstract 

In this article, we study the interactions of stable, hadronizing new states, arising in certain 
extensions of the Standard Model. A simple model, originally intended for stable gluino hadrons, 
is developed to describe the nuclear interactions of hadrons containing any new colour triplet 
or octet stable parton. Hadron mass spectra, nuclear scattering cross sections and interaction 
processes are discussed. Furthermore, an implementation of the interactions of heavy hadrons 
in GEANT 3 is presented, signatures are studied, and a few remarks about possible detection 
with the ATLAS experiment are given. 
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1 Introduction 



Among the more plausible scenarios of physics beyond the Standard Model is supersymmetry 
(SUSY) pQ. In conventional SUSY models, the lightest supersymmetric particle (LSP) is neutral 
and colourless. However, models exist in which the LSP is coloured. Models with a gluino LSP 
are reviewed in Ref. [2], and include gauge mediated supersymmetry breaking models and string 
motivated supersymmetric models. Stable squarks are cosmologically disfavoured because of their 
nonzero electric charge [3], but could be sufficiently long-lived to behave like an effectively stable 
LSP in a detector experiment. A coloured LSP would hadronize into heavy (charged and neutral) 
bound states. These bound states (for example gg, gqq, gqqq, qq, qqq) are generically called R- 
hadrons, where the "R" refers to the fact that they carry one unit of R-parity In this article, we 
shall focus on the fact that such a hadronized LSP will have measurable interactions in a detector, 
in contrast to conventional SUSY studies, in which the LSP is a weakly-interacting particle, i.e. a 
neutralino or a sneutrino. 

In addition to supersymmetry, other extensions of the Standard Model have been proposed, which 
predict the existence of new heavy hadrons, either due to the presence of a new conserved quantum 
number, or because the decays are suppressed by kinematics or couplings. For example, theories 
with leptoquarks predict stable hadronized states, if the Yukawa coupling between the leptoquark 
and its decay products is so small that the leptoquark is long-lived |5J- Another example is theories 
with universal extra dimensions, where exact momentum conservation in all dimensions leads to 
stable Kaluza-Klein excitations of e.g. quarks and gluons, which would form stable hadronizing 
states 

Searches for stable massive particles are summarized e.g. in Ref. [Jj, which includes different search 
methods. Negative results of searches in ordinary matter and cosmic rays imply that the relic 
density of such particles must be too small for detection. In Ref. [2J it is argued that an LSP 
gluino can indeed have a very small relic density. Alternatively, the cosmological bounds on stable 
massive particles can be evaded if the hadron decays slowly. In this paper, we refer to particles as 
stable as long as they are stable from a detector point of view. Concerning search methods using 
accelerators, these are limited by the center of mass energy of the accelerator and are concentrated 
on charged heavy particle searches ,8... A few searches for neutral hadrons, in particular R-hadrons, 
exist |10(lll|. From accelerator searches it is clear that, if heavy hadrons with masses less than the 
order of a hundred GeV had existed, they would have been detected. In this paper, particles in 
which we are interested, and which we classify as heavy, are particles with masses > 100 GeV. 

The energies of heavy hadrons produced at a p — p collider like the LHC are typically about two to 
three times their own mass, i.e. they may be relativistic but their mass is still far from negligible. 
On their trajectory through the detector, they slow down even further. 

Although electromagnetic interactions of heavy states are well understood, little information is 
available about nuclear interactions of heavy hadrons, and the resulting energy deposits in a detector 
calorimeter. Models to describe the nuclear interactions of gluino bound states in particular have 
previously been discussed in Refs. [21 1121 ITS] . Here, a new model for the interactions of gluino 
hadrons is presented, which is more generally applicable to any kind of heavy hadrons, and which 
provides a convenient basis for future refinements. 

The organization of this paper is as follows. First, electromagnetic interactions of charged heavy 
hadrons will be discussed in Section |2J Next, in Section 3, we propose a model describing nuclear 
interactions of heavy hadrons. Aspects such as hadron mass spectra and nuclear scattering cross- 
sections will be discussed. Section 4 is devoted to the GEANT 3 simulation of R-hadronic nuclear 
interactions. In Section 03 results of the simulation will be presented, together with studies of 
typical detector signals. In Section |SJ we briefly discuss manifestations of R-hadrons in the ATLAS 
experiment, and finally in Sectional we summarize and discuss future studies. 
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2 Electromagnetic interactions of heavy charged hadrons 



A charged heavy particle suffers both continuous ionization losses as well as repeated Coulomb 
scatterings [§]. Continuous ionization losses of heavy particles are known to be described by the 
Bethe-Bloch equation. Since electromagnetic losses are proportional to 1//3 2 , the losses for a particle 
moving with (3 <Cl are considerable. The ultra-relativistic rise in electromagnetic energy losses is 
not relevant here, as can be seen from Fig El Repeated Coulomb scatterings change the particle 
trajectory. The deflection angle, being proportional to l//3p, is small, the small velocity being 
compensated for by the large momentum with which a heavy particle is normally produced. 

3 Nuclear interactions of heavy hadrons 

Both charged and neutral heavy hadrons lose energy through scattering off nuclei. In the following, 
a simple and general framework is presented for simulating nuclear interactions of heavy hadrons, 
independent of the new physics model in which the hadron arises. Exotic colour triplet states will 
be denoted by C3, colour antitriplet states by Cg and colour octet states by Cg. When referring to 
the mass of a parton, we always mean the constituent mass. An interpretation in terms of a pole 
mass or other mass definitions is clearly irrelevant for the qualitative considerations in this paper. 
For simplicity, we take into account only states containing u and d quarks, with constituent masses 
assumed to be m q = 0.3 GeV. 

3.1 The role of the heavy object 

Any stable heavy coloured exotic particle hadronizes and forms a colour singlet, for example C3C;, 
Csqq, Cgq, Csqq, Csqqq, Csg, etc. The probability that the parton Q of colour state i will interact 
perturbatively with the quarks in the target nucleon is small, since such interactions are suppressed 
by the squared inverse mass of the parton. As a consequence, the heavy hadron can be seen as 
consisting of an essentially non-interacting heavy state Cj, accompanied by a coloured hadronic 
cloud of light constituents, responsible for the interaction. The effective interaction energy of the 
hadron is therefore small, as can be seen by considering a Csqq state, e.g. with a total energy 
E=450 GeV and a mass m of the Cs parton of 300 GeV, i.e. with a Lorentz factor of 7=1.5. 
Although the kinetic energy of the hadron is 150 GeV, the kinetic energy of the interacting qq 
system is only (7 — l)m. qq ~ 0.3 GeV, (if the quark system consists of up and down quarks). Thus, 
the energy scales relevant for heavy hadron scattering processes off nucleons are low! It is therefore 
most likely that interaction processes are mediated by Reggeons (Fig. ^1) and not by by Pomerons 
(Fig. In conclusion, it becomes apparent that the presence of the heavy parton Cj has two 
basic consequences. 

• It acts as a reservoir of kinetic energy. After an interaction, where the light interacting system 
loses kinetic energy, new kinetic energy is transferred to it from the co-moving Cj-parton. 

• The parton Cj forces the quark system to be in a certain colour-state, in order for the system 
as a whole to form a colour-singlet state. This influences the mass spectrum of the hadrons, 
as will be explored in detail in Section f3. 21 

The following issues are not influenced by the presence of the heavy parton: 

• The interactions of the light constituents are of the same character as for an ordinary hadron. 
As stressed in the beginning of this section, the heavy state Cj acts as spectator in an 
interaction. 
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Figure 1: R-hadron scattering of a proton, (a) Reggeon mediated elastic scattering, (b) Pomeron 
mediated elastic scattering, (c) the formation of nuclear resonances 



• The overall size of the hadron as a whole is probably not strongly influenced. This can 
be understood by recalling that the wave-function associated with a particle with mass M 
roughly scales as 1/M 2 . The influence of the heavy object on the size of the total hadron 
is therefore minimal, and the light quarks dominate completely. 



3.2 Hadron mass spectrum 

In order to establish the signatures of heavy hadrons in a detector, the mass spectrum of the 
states must be understood. For example, if the Csuuu state were to be the lightest baryonic 
state, then decays into it would be kinematically favourable. This would be important, since the 
electromagnetic energy losses would then be increased significantly (if the Cs parton does not carry 
negative charge). A completely different signature would result if a neutral Csudd state would be 
the lightest. 

The mass of the lowest-lying hadronic states, with no radial excitation or orbital angular momen- 
tum, can be approximated by (Ref. [Tl| ) 

E, ( F i • F j)(Si • Sj) 
rrii-ky =! — (1) 
mirrij 

i i^j ■> 

in which the summation is over all partons i contained in the hadron, mj is the parton constituent 
mass, Fj is the SU(3) colour matrix for parton i, Si is the SU(2) spin matrix for parton i, and k is 
a constant with dimension (mass) 3 . The second term, responsible for the mass splitting, depends 
on the colour-state and spin-state of the hadron. In the derivation of the hadron mass spectrum of 
heavy hadrons, terms involving the heavy parton in the denominator can be neglected. Expressions 
for (Fi • F2) and for (Si • S2) can easily be derived for a qq and qq state [TIJ. A qq state can either 
form a colour antitriplet or a sextet configuration (3 (g> 3 = 3 © 6), while a qq state can form a 
colour singlet or an octet state (3<8> 3 = 1©8). The eigenvalues of the operator (Fi -F^) are —4/3, 
—2/3, 1/3 and 1/6 when the state is in a colour singlet, antitriplet, sextet and octet configuration, 
respectively. The eigenvalues of (Si • S2) are —3/4 and +1/4 for a qq or qq state with total spin 
zero and spin one, respectively. Eq. ^ leads to the familiar mass patterns of the known hadrons. 
For example, the tt — p, K — K*, D — D* and B — B* systems follow Eq. n reasona bly well, with 
k ps 0.043 GeV 3 . Below, the mass spectra of heavy meson and baryon states will be discussed using 
Eq. ^ With more refined methods, it would be possible to calculate the spectrum of all possible 
heavy mesons and baryons precisely. Here, we are mainly interested in the order of magnitude of 
the mass splittings and not in fine details. 
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3.2.1 Meson mass spectrum 



Possible heavy mesons are Cgq, C^q and Csqq states. For Cgq and C%q states, no significant mass 
splitting is expected to occur, since the second term of Eq. ^ is negligible (certainly significantly 
smaller than that of the B — B* system, which is only 46 MeV). 

For a Csqq state, the mass spectrum is given by: 

M C8qq « M Cs + 0.3 + 0.3 - 0.043 x fj^^ s qq =0 
ps M Cg + 0.66 

M C8qq ~ M C8 +0.3 + 0.3 - 0.043 x s qq =l 
ps M Ca + 0.58 

There are two noticeable aspects. First, the mass hierarchy is reversed, as compared to that of the 
7r — p mass splitting: the spin-zero state is the heaviest. Second, contrary to the ir — p case, the 
mass splitting between mesons with different spin is much smaller. Consequently, mass splittings 
for heavy mesons may safely be neglected. 



3.2.2 Baryon mass spectrum 

Possible heavy baryons are C3qq and Cgqq, Csqqq and Csqqq states. The mass spectra are obtained 
by a similar calculation as for the mesonic case, with k ps 0.026 GeV 3 , as derived from the ordinary 
baryon sector. For C3qq states (and, by symmetry for the Cgqq baryons) we obtain the mass 
spectrum 

C_2 x _3n 

M C3qq « Mc 3 + 0.3 + 0.3 - 0.026 x q 3 x q 3 s qq=° 
ps M Ca + 0.46 

(-- X +i) 

M C3qq « M C3 + 0.3 + 0.3 - 0.026 x^|—^ s qq =l 
ps M Ca + 0.65 

At this stage we must recall that the total wavefunction associated with a quark system can be 
decomposed in flavour x spin x colour, which has to be anti-symmetric for the qq system here. 
Taking into account that the colour wavefunction, associated with the two quarks, is antisymmetric 
(it is in an antitriplet configuration), it is seen that the two quarks in a C3UU or Cgdd state can 
only be in a symmetric spin-configuration, i.e. s qq = 1. This implies that it will kinematically be 
favourable for the other baryons to decay into the lighter Cgud state with s qq = rather than in 
the heavier C3UU, Cgud or Csdd states with s qq = 1. 

Deriving the mass spectrum of Csqqq is considerably more complicated. Fortunately, however, we 
expect the mass splitting to be small. This is based on the following line of arguments. 

Recall that the total wavefunction associated with the three quarks in a Csqqq state must be 
anti-symmetric. For the Csuuu or Csddd states, the three quarks are in a symmetric flavour 
configuration, and thus the spin x colour wavefunction should be anti-symmetric. The three quarks 
being in a colour octet configuration implies that s qqq =l/2 14 j. Thus, the qq states involved in 
Eq. [I] have spin s qq = 1 or s qq — and the associated wavefunctions are mixed symmetric and 
anti-symmetric. The colour configuration of the qq states is also mixed: either colour antitriplet 
or sextet. 

For the Csuud and Csudd states, the situation is slightly different. The wavefunction associated 
with flavour can either be totally symmetric or mixed symmetric and anti-symmetric. In the 
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former case, the situation is exactly the same as for the Cguuu and Cgddd states. In the latter 
case, the spin x colour wavefunction associated to the three quarks must also be mixed. The three 
quarks being in a colour-octet configuration, such mixed states can be obtained with s qqq =l/2 
or Sqq q = 3/2 ^1] . The spin doublet case, s qqq =l/2, implies a mixed colour configuration and a 
mixed spin configuration for the three quarks together. For the quark pairs involved in Eq. this 
means that they have a mixed s qq = 1 or s qq = configuration, and a mixed colour antitriplet 
and sextet configuration. On the other hand, the spin quartet case, s qqq = 3/2, implies a mixed 
colour configuration and a symmetric spin configuration for the three quarks. For the quark pairs 
involved in Eq. ^ this means that they have a symmetric s qq = 1 spin configuration, and again 
their colour configuration is a mixture of colour antitriplet and sextet. 

Summarizing, Cguuu or Cgddd states only arise if s qqq =l/2, while Cguud and Cgudd states arise 
if s qqq =I/2 or 3/2. (Note that this is the opposite for ordinary baryons!) We saw above that 
both in the s qqq = 1/2 case and in the s qqq = 3/2 case the wavefunctions of the qq combinations 
involved in Eq. ^are always strongly mixed: the colour wavefunction is mixed, and either the spin, 
or the flavour wavefunctions, or both, are mixed for the qq combinations. The colour antitriplet 
and sextet contributions enter with opposite sign, as do the s qq = 1 and s qq = terms. Without 
having done the exact calculation, it is therefore certain that this feature will result in a partial 
cancellation of the mass splitting expression. Thus, the splitting is not expected to be more than 
the order of a hundred MeV. Exactly the same applies for Cgqqq states. In conclusion, the mass 
splittings between the different Csqqq and Csqqq spin states are not expected to be large and may 
therefore be neglected to first approximation. 

3.3 Resonances 

The formation of resonances is closely connected with the discussion above. In the presence of a 
colour state Q, new resonances may arise, as illustrated in Fig. where Q = Cg is a gluino. We 
do not have any knowledge about the detailed mechanism of resonance formation in heavy hadrons, 
but a few remarks can be made. Concentrating on Cg hadrons, if indeed mass splittings between 
the different hadrons are not larger than the order of a hundred MeV, resonances will not play a 
significant role. The minimum center-of-mass energy for a scattering of for example a Cgqq state 
with a nucleon with mass m n is mc 8 qq + ni n mc a + mc s + 1.5 GeV, an energy which is 

above the masses of the most significant resonances. Thus, even if resonances exist, they are not 
likely to play an important role, but they would be smoothly merged with the continuum, as for 
ir — p scattering above the A resonance region. Because of this, and lacking a detailed description 
of these resonances, we will not take them into account explicitly. 

3.4 Black disk approximation 

Predicting the total cross section of a heavy hadron scattering off a nucleon is non-trivial. Here, 
a series of arguments will be presented, that will be used as guidelines in constructing an effective 
model. At high center-of-mass energies, cross sections may be approximated by the geometrical 
cross section. At low energy, the cross section is hard to estimate; in fact, even for normal hadrons, 
low energy scattering cross sections are poorly understood. A vanishing cross section at low energy, 
as exhibited by pion-nucleon scattering, is not probable. The reason why this cross section vanishes 
is that s-wave scattering is forbidden, a unique feature for pion-nucleon scattering [T^]. In a Cgqq 
hadron, even though containing u and d quarks, the quark system is in a colour-octet state, and 
is thus no ordinary pion. In particular, a colour octet pion is not light, as shown in the previous 
section, and can consequently not be treated as a Goldstone boson; hence s-wave scattering should 
be possible. On the other hand, a rise in the cross section in for example the nuclear scattering of a 
Csqqq state at low energies, as is the case for proton-neutron scattering, is likewise also improbable. 
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The reason why this cross section is so large is again a unique one and is connected to the existence 
of the deuteron |15j . Expecting neither a rise, nor a vanishing behaviour, the pragmatic solution 
is to treat a heavy hadron simply as a black disk and to use the geometrical cross section at all 
scattering center-of-mass energies. 

3.5 The size of the total cross sections 

The size of the heavy hadron being roughly the same as the size of the accompanying hadron system, 
the total cross section for nucleon scattering can be approximated by the asymptotic values for the 
cross sections for normal hadrons scattering off nucleons. For example, in the case of a Csqq state, 
in which the qq system contains only u and d quarks, one can use the asymptotic value of the 
total cross section for pion- nucleon scattering, while in case one s quark is present, one can use 
asymptotic values for kaon-nucleon scattering, etc. A simple rule would be 12 mbarn for every u or 
d quark, and 6 mbarn for every s quark scattering off a nucleus. Cross sections for any Q-hadron 
can then be calculated. As said before, we take into account in this article only u and d quarks. 
Although it is known that a significant amount of s quarks should in principle be present in the 
gluino R-hadrons (ca 15%), only a small error on the cross section is made, the difference only 
being 6 mbarn. The amount of hadrons with two s quarks or heavier is negligible. For a gluino 
R-meson and R-baryon, the total cross sections are 24 mbarn and 3/2 x 24=36 mbarn, respectively. 
A gluino-gluon state can be assumed to have the same cross section as a gluino R-meson, since 
the geometrical cross section is approximated by the high energy hadron cross section, where gluon 
exchange would dominate. The gluon-gluon coupling is a factor 9/4 larger than the quark-gluon 
coupling, but a meson has two quarks, resulting in a cross section of a gluino-gluon state which is 
(9/4)/(l+l)~l times the cross section for a gluino R-meson. 

3.6 The size of the partial cross sections 

Besides the total cross section, which determines the interaction length, the partial cross sections 
for 2 — ► 2 and 2^3 processes must be known in order to estimate the energy loss per nuclear 
interaction. An estimate of the number of final-state particles can be made by studying the energy 
Q available for the production of kinetic energy and potentially extra particles. In the case of a 
2 — * N scattering, where an incoming particle with mass mi scatters on a particle at rest with 
mass m,2, we define the Q-value by 



where rrij is the mass of final-state particle i. The value of s/s is mf + m\ + 2E / f b m2, where E^ b is 
the energy of the incoming particle. If the Q-value exceeds the mass of a pion, then an extra pion 
may be produced. To determine the relevant order of magnitude of Q in the scattering of heavy 
objects, consider an elastic scattering of a heavy particle with mass mi, moving with a Lorentz 
factor 7, off a light particle at rest with mass m^- It is easily seen that ^/s ~ mi + 77772 and 
Q (7 — l)m.2, the latter being obviously small if the incoming heavy particle has a small value 
of 7. Fig. |2ji displays the Q-value of a process R + n — ► R + n, in which a heavy gluino hadron R 
with mass M scatters off a nucleon n with mass m n at rest. Fig. \^ shows that for heavy hadrons, 
the Q-value is much smaller in comparison with light hadrons of the same momentum, implying 
that the number of final-state particles is also much smaller. For the range of momenta of particles 
produced at a collider like LHC, the Q values are so small that the final-state multiplicity is rarely 
expected to exceed three. Even though the Q-values may exceed the pion mass, observation of 
for example pion-proton scattering shows that three final-state particles are produced only above 
yfs ~ 1.5 GeV, where Q « 0.4 GeV 3> m n . The relative amounts of the 2 — > 2 and 2^3 processes 




(2) 
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Figure 2: (a) The value of Q = yfs — m n — ran, f° r hadrons with mass M scattering with lab 
momentum on a nucleon at rest, (b) The value of the cross section for gluino R-meson-proton 
scattering (horizontal lines) in comparison with pion-proton scattering, as well as the determination 
of Qq and its uncertainty (curved dashed lines connecting the horizontal lines). 



occurring for heavy hadrons at high energy can be derived from high energy pion-proton scattering 
data, where the cross section for 2^2 and 2^3 processes is known to be 15% and 85% of the 
total cross section, respectively. In Fig. [2j3, such cross sections are displayed for gluino R-meson 
nucleon scattering, together with the cross sections for pion-proton scattering. At low momentum, 
the amount of 2 — * 3 processes is then obviously overestimated, and sometimes 2 — ► 3 processes 
are even kinematically impossible. To solve this, a phase space function describing the relative 
amounts of 2 — ► 2 and 2 — ► 3 processes is introduced. 



3.7 Phase space considerations 

In order to predict accurately the energy loss of heavy hadrons in matter, the relative amount of 
2 — * 2 and 2 — ► 3 processes must be estimated. A constant matrix element is assumed for both 
processes, so that the relative yield of the three-body final states is initially determined solely by 
a phase space factor. A function describing the available phase space is expected to be a function 
of Q, defined in Eq. EJ and a reasonable ansatz is given by 

F ^ = d$ 2 (Q) d$7(Qp ( 3 ) 
d-MQo) d$ 3 (Qo) 

where d<3? n is the available n-body phase space. The fact that the two- and three-body phase spaces 
do not have the same dimension (compensated for by the here neglected matrix elements) forces 
us to normalize phase space to that available at a certain value Qq. For Q < ra^, the phase space 
for three-body scattering is vanishing, while for Q —* oo, it is entirely open. Simplified expressions 
for d$2 and d$3 can be derived under the assumption that the mass of the heavy hadron is much 
larger than the nucleon mass. If one of the final-state particles in the three-body processes is a 
pion, a general expression describing the relative amount of three-body phase space as function of 
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(a) (b) 

Figure 3: (a) The phase space function F(Q) drawn for a few values of Qq and (b) as a function of 
the lab momentum Pi a b of a heavy hadron with mass M scattering off a nucleon at rest, and with 
Qo = 1.1 GeV. 



Q, if Q > m^, may be derived: 

F(Q) = ^ {l ^3 )V2 . (4) 

W( 1+ ^)© 3/2 

From the optical theory for light hadrons, it is known that elastic scattering accounts for roughly 
15% of the total cross section at high momenta, as is already mentioned in the previous section. 
Taking this into account, a function f{Q) representing the relative probability for 2^3 scattering 
is 

f(Q) = 0.85F(Q) (5) 



An appropriate value for Qq can be determined experimentally, e.g. by finding the value where 
F(Q) is exactly 1/2, thus where f(Q) = 0.425. Fig. illustrates the way in which we determine 
empirically the value for Qq and its uncertainty. Three smooth dashed curves are drawn in between 
the horizontal lines corresponding to the 2^2 and the total cross sections. The central curve is the 
best estimate of the division between the elastic and inelastic cross section behaviour neglecting 
resonances, while the other two curves represent the uncertainty. The point in which F(Q) = 1/2 
corresponds to a purely elastic cross section of (0.15 + 0.425) x 24 mb = 14 mb, which occurs at 
yfs 1.7 GeV in the central curve, where Q = 1.7 — m p — m n = 0.6 GeV/c 2 . The lower and upper 
limit for this point are Q = 1.4 — m p — m n = 0.4 GeV/c 2 and Q = 2.1 — m p — m n = 1.1 GeV/c 2 . 
Substitution of Q in Eq. |1] gives Q = 1.1, with upper and lower limits of 0.6 and 1.7 GeV/c 2 , 
respectively. This value roughly fits with the values obtained for p — p and p — i^ + -scattering 
processes. The resulting phase space function F is drawn in Fig. as function of Q, as well as 
the upper and lower limits. In Fig. |3t> , the phase space function with Qq = 1.1 GeV/c 2 is drawn 
as function of the lab momentum of the incoming hadron for different mass values. 
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3.8 Nuclear scattering processes 



As discussed in Section |3~T1 the presence of the gluino changes the colour-state of the hadron cloud, 
but not the character of the interactions of this cloud. Since scattering processes are low-energetic, 
all scattering processes for mesons and baryons can simply be derived by exchanging quarks. Below, 
interactions of hadronized colour octet and colour triplet states will be discussed. 

3.8.1 Interactions of Cgqq 

Interactions of Csqq states include the following processes. The 2 — > 2 processes are purely elastic 
scattering (e.g. Csdd + uud — ► Csdd + uud), charge exchange (e.g. Csdd + uud — ► Csud + udd) 
and baryon exchange (e.g. Csdd + uud — * Csudd + ud) while the 2^3 processes include normal 
inelastic scattering (e.g. Csdd + uud — ► Csud + udd + dd) and inelastic scattering with baryon 
exchange (e.g. Csdd + uud — ► Csuud + ud + du). 

It should be noted that the processes with baryon exchange are kinematically favoured due to the 
fact that a final-state pion is so light. In such processes, typically an extra 2m^ — m n ~ 500 MeV 
of kinetic energy would be liberated. However, these processes could be dynamically suppressed 
because two quarks must be exchanged. 

3.8.2 Interactions of Cgqqq states. 

Interaction processes include 2 — ► 2 processes like purely elastic scattering (e.g. Csuud + uud — > 
Csuud + uud) and charge exchange (e.g. Csuud + udd — > Csudd + uud) and 2 — ► 3 processes 
like Csuud + udd — > Csudd + uud + dd. No baryon exchange is possible, since the probability 
that the Csqqq state interacts with a pion in the nucleus is negligible, and besides this process 
would kinematically be strongly disfavoured. This has the important implication that mesons get 
converted into baryons during repeated interactions, but not vice versa. 

3.8.3 Interactions of Csqqq states. 

Only a very small amount of Csqqq states will arise in the hadronization process. Interactions differ 
from those of Csqqq states, since scattering takes place on protons or neutrons, containing quarks 
rather than antiquarks. A Csqqq state would thus interact dominantly by baryon annihilation, for 
example Csuud + uud — > Csud + ud. This process would kinematically be favourable. 

3.8.4 Interactions of Cgg states 

A gluon is able to convert into a uu or dd state, and as such, a Csg state probably interacts like 
(and mixes with) Csuu or Csdd states. The mass of the active system, the gluon, is usually taken 
to be 0.7 GeV approximately the same as the constituent mass of two first generation quarks. 
Thus, the possible interaction processes are expected to be similar to those for Csuu or Csdd 
states. Besides the fact that the interaction processes are similar, the cross section is expected to 
be roughly the same, as has already been explained in Section 1331 allowing us to treat a Csg state 
like a neutral Csqq state. As for Csqq states, the Csg states will eventually convert into baryons. 

3.8.5 Interactions of C3C1 and Cgq states 

Interactions of a C3CI state include processes with quark-antiquark annihilation. Thus, there are 
2 — > 2 processes such as elastic scattering (e.g. C3U + uud — ► C3U + uud), charge exchange (e.g 
C3C1 + udd — > C3U + uud), or baryon exchange (C3U + udd — > Csud + iid). For 2^3 processes, 
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these include normal inelastic processes, like C3U + udd — > Csd + udd + ud, as well as processes with 
baryon exchange, like C3U + udd — > Csud + ud + uu. As above, baryon exchange is kinematically 
favoured by the possibility of having a light pion in the final state. 

Interactions of a Cgq state include similar processes but without quark-antiquark annihilation and 
without baryon exchange. Possible processes in that case are 2 — > 2 processes like elastic scattering 
(e.g. Cgu + uud — > Cgu + uud), charge exchange (e.g. C§d + udd — ► C§u + uud), and 2^3 
processes like Cgu + udd — > Cgd + uud + dd. 

3.8.6 Interactions of C3qq and C§qq states 

Interactions of C3qq states are processes without quark annihilation and include 2^2 processes 
with purely elastic scattering (e.g. C3UU + uud — > C3UU + uud) and charge exchange (e.g. C3UU + 
udd — > C3ud + uud), and 2 — > 3 processes like e.g. C3UU + udd — ► C^ud + udd + ud. 

On the other hand, Cgqq states may interact by quark annihilation and thereby baryon annihilation. 
Processes include 2 — ► 2 processes like purely elastic scattering (e.g. Cguu + uud — ► Cguu + uud), 
charge exchange (e.g. Cgiiu + uud — > Cgud + udd), and baryon annihilation (e.g. Cgiiu + uud — > 
Cgu + ud) and 2 — > 3 processes like e.g. Cguu + udd — ► Cgud + udd + ud. Baryon annihilation 
would kinematically be favoured. 

3.9 Relative probabilities of scattering processes 

In general, an enormous number of scattering processes is possible, over 140 when summed over 
all gluino R-hadrons. To know which processes take place, the target (neutron or proton) must 
be known, as well as the relative coupling of all processes. The latter requires the calculation of 
the Clebsch-Gordan coefficients of isospin-related processes, plus an evaluation of all additional 
dynamical effects for all processes. This however has not been done. Since the masses of the 
different lowest-lying Cs mesons are degenerate, and similarly for the baryons, as discussed in 
Section l3~2l the phenomenology is not expected to be affected by this assumption. Therefore, after 
the target and the class of interaction processes, i.e. 2 — > 2 or 2 — > 3 processes, is determined, equal 
relative couplings for all processes inside this class are assumed. 

3.10 The nuclear cascade 

The incoming hadron interacts with a nucleon inside a nucleus. Issues like the Fermi motion of the 
nucleons inside the nucleus, binding energy of the nucleus, evaporation energy of a nucleus, and 
instability of the nucleus after the interaction play an important role. It depends on the energy 
of the incident hadron whether the interaction causes a large nuclear cascade or not. To estimate 
this, the wavelength associated with the incident hadron, A = can be compared with the size 
of the nucleus. Since the interacting system of a heavy hadron is low-energetic, the associated 
wavelengths are of the order of the size of the nucleus, and the development of a nuclear cascade is 
probable. 

3.11 Nuclear energy losses 

As argued in Section ESI the total energy loss in a collision of the heavy hadron is small. This small 
energy transfer is shared between the kinetic energy of the kicked-out nucleon or pion, its binding 
energy in the nucleus, the production of extra final-state particles, and nuclear degrees of freedom 
inside the remnant. 
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Figure 4: The tracking steps in GHEISHA. 

4 GEANT 3 simulation of R— hadronic interactions 

The simulation of processes which accompany the propagation of an R-hadron through the detector 
material is performed within the framework of GEANT 3. For the specific case of describing R- 
hadronic interactions in matter, the GHEISHA package is used ^7j- Dedicated routines to describe 
interactions of gluino R-hadrons are developed |18j . Tracking of R-hadrons through a volume is 
done without modification. For the generation of R-hadronic nuclear interactions, the existing 
pion routines are used as starting point, following the idea of Refs. ^Hl and JU]. In Fig. @J the 
GHEISHA tracking steps and the new routines are displayed. Three main issues characterize the 
generation of hadronic interactions in GHEISHA, which are simulated as follows. 

1. The evaluation of the cross sections, needed to calculate the mean free path. Cross sections 
are calculated according to the arguments in Section 13.51 

2. The selection of the interaction and subsequent sampling of the final-state multiplicities. 
The target and selection of the interaction process is done with the help of seven new shower 
routines for gluino R ++ ,R + ,R° and R~ baryons and for R + ,R° and R~ mesons. Due to 
the tiny amount of Csqqq states no shower routine routines devoted to this states have been 
developed. The interaction processes are chosen according to the prescription in Section G3 

3. The generation of the final-state particles and their kinematics. The pr-values for the final- 
state particles are tabulated from experimental data. For R-hadrons, this implies that t values 
are selected as in the pion routines following da / dt=e~ bt , where b is an empirical function of 
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Figure 5: The fraction of R-mesons converted into R-baryons as function of the traveling length, 
starting with a neutral R-meson. 

the lab momentum of the incoming hadron. A rescaling of momenta to the momenta of the 
active quark (or gluon) system is applied for R-hadrons. 

The fact that the hadron scatters on a nucleus and not on a free proton is taken into consideration in 
the three issues mentioned above. Apart from a rescaling of the evaporation energy of the nucleus, 
being a function of the kinetic energy of the incoming hadron, to the kinetic energy of the active 
incoming quark (or gluon) system, the issues are simulated exactly as for pions. The code, as well 
as more detailed information about the simulation, is available in Ref. 18]. 



5 Results 

5.1 Meson-baryon conversion 

As explained in detail above, an R-meson is able to convert into an R-baryon by scattering off a 
nucleon, but not vice-versa. The amount of converted mesons as a function of the traveling length 
in a piece of iron is shown in Fig. OD Note that an R-baryon has a larger cross section for nuclear 
interactions than an R-meson. The conversion therefore increases the subsequent energy losses in 
a calorimeter. 



5.2 Kinematics 

We cross-check the values for the scattering angle in the center of mass system, cos 6*, and for 
the transverse momentum pt of the scattered R-hadron. Results are shown in Fig. At small 
momenta, the 2^2 scattering processes are isotropic, while at large momentum, a clear forward 
peak is seen, explained by the fact that the pt in nuclear scatterings is limited to approximately 1 
GeV. 
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5.3 Losses per collision 

To check whether non-trivial nuclear effects are simulated correctly, the losses per nuclear collision 
have been evaluated in iron, carbon and liquid hydrogen. In Fig. the losses per collision are 
displayed for 2^2 and 2^3 scatterings in the three materials for different values of the R- 
hadron mass. From the three figures obtained for e.g. iron, it can be seen that the loss per collision 
decreases with increasing R-hadron mass for fixed kinetic energy. This is logically explained by the 
fact that the energy available for the production of new particles in a scattering becomes smaller. 
From these plots, we also note that the energy loss per collision is generally much larger for a 
2 — ► 3 processes than for 2 — ► 2 processes, which follows from kinematics. If we compare different 
target materials, the loss per collision increases for heavy elements. This can be explained by the 
fact that the R-hadron is slow, and can scatter several times with the different nucleons inside one 
nucleus. For hydrogen, it should be noted that the energy loss for 2 — ► 3 processes vanishes at small 
momenta, a feature which is not seen for heavier elements, where this effect is washed out by the 
nuclear cascade. This is due to a technical problem in the simulation package GHEISHA, which 
is entirely based on parameterizations for pions. Repairing this 'bug', and making the package 
optimal for heavy hadrons, would require a considerable rewrite of the code. This is beyond the 
scope of this work and, besides, 2^2 scattering processes dominate completely at low energies. 
Energy losses for an R-hadron which converts from being a meson into baryon are slightly smaller 
when considered from the R-hadron perspective, but the total energy deposited in the detector is 
larger due to the larger kinetic energy released in the process. 

5.4 Nuclear energy losses in different materials 

We compare the energy loss to purely nuclear interactions (i.e. no ionization losses) of an R-hadron 
in different materials, when the total amount of nucleons traversed is the same. Total nuclear 
losses of R-baryons in 10 cm iron, 36 cm carbon and 11m liquid hydrogen are shown in Fig. |H1 
It turns out that approximately the same amount of total energy is lost. The average number of 
nuclear interactions is however not the same. As can be easily derived from the cross section on 
free nucleons, the interaction lengths in iron, carbon and liquid hydrogen of an R-baryon are 14 
cm, 32 cm and 610 cm, respectively. In 10 cm iron, 36 cm carbon, and 11 m liquid hydrogen ca 
0.7, 1.1 and 1.8 collisions take place. However, the amount of energy lost in a collision increases 
with increasing atomic number (see Fig. EJ), and the overall result is that approximately the same 
amounts of energy are lost. 
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Figure 6: At the left, distributions for cos#* in 2— >2 scattering for different values of the R-hadron 
initial momentum (left). At the right, scattering plots for the generated px value of the R-hadron 
in 2^2 scattering for three different values of the R-hadron mass: 100, 300 and 600 GeV/c 2 . 

14 




Momentum (GeV/c) Momentum (GeV/c) Momentum (GeV/c) 



> 
o 



Iron 



















200 400 600 800 
Kinetic energy (GeV) 



> 

a 



Carbon 



200 400 600 800 
Kinetic energy (GeV) 



> 

a 




200 400 600 800 

Kinetic energy (GeV) 




200 400 600 800 
Kinetic energy (GeV) 



r-, 2 
> 

a 

c 

.2 1.5- 



Carbon 



200 400 600 800 
Kinetic energy (GeV) 



> 

a 



0.2 



Hydrogen 



:- 



200 400 600 800 
Kinetic energy (GeV) 




200 400 600 800 
Kinetic energy (GeV) 



1.5 

> 



Carbon 











200 400 600 800 
Kinetic energy (GeV) 



> 

a 



0.6 - 



Hydrogen 











200 400 600 800 
Kinetic energy (GeV) 



Figure 7: Profile plots representing the mean value of the average energy loss per collision of an 
R-baryon for three different values of the R-hadron mass: 100 (three upper Figures), 300 (three 
middle figures) and 600 GeV/c 2 (three lower figures). The lines represent the mean value and 
its error. The lower and upper profiles in each plot correspond to the energy loss in a 2 — * 2 
scattering and a 2 —* 3 scattering, respectively, while the middle profiles represent the average loss 
per collision. 



5.5 Total energy losses 



The total energy loss in 1 m iron, i.e. the sum of nuclear energy losses and ionization losses, is 
determined for different mass values of the R-hadron, and results are displayed in Fig. The 
ionization losses for a singly charged R-hadron traversing 1 m iron without suffering any nuclear 
interactions are shown as well. Nuclear losses are seen to be completely dominating at high energies. 
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Figure 8: Comparison of the total nuclear energy losses in 10 cm iron, in 36 cm carbon and in 11 
m liquid hydrogen, corresponding to the same amount of nucleons, for an R-hadron mass of 100 
GeV (upper curves), 300 GeV (middle curves) and 600 GeV (lower curves). 



5.6 Nuclear energy losses versus ionization losses 

A direct comparison between the energy losses of an R-hadron suffering only nuclear collisions, and 
those of an R-hadron suffering ionization losses as well as nuclear collisions, is shown in Fig. 1101 
At a kinetic energy of 300 GeV, ionization losses are responsible for a shift of the mean value of the 
energy loss of 0.9 GeV, slightly less than the value displayed in the ionization figure for M = 100 
GeV in Fig. £0 This is obviously due to the fact that the R-hadron may change charge and hence 
may be neutral during part of the trajectory. At a kinetic energy of 20 GeV, the shift is larger, 1.9 
GeV for this particular R-hadron mass. The conclusion is that for low-energetic heavy hadrons, 
ionization losses are comparable to nuclear losses, while at high kinetic energies, nuclear losses are 
completely dominating, as expected. 



5.7 Dependence of energy losses on the phase space function 

To evaluate the dependence on the parameter Qq, appearing in the phase space function F(Q), 
it is varied from 0.6 GeV/c 2 to 1.7 GeV/c 2 . From Fig. El we note that the available phase space 
for 2 — > 3 scattering decreases slightly with increasing value of Qo, and therefore the energy loss 
is expected to decrease for increasing Qq. For two values of the R-hadron mass, the total energy 
loss in 1 meter iron is studied for different values of Qo, and results are shown in Tabled The 
dependence on Qo of the total energy loss is small, and indeed losses are smaller for increasing 
Qo-value. 



5.8 Penetration depth 

In order to estimate the penetration depth of R-hadrons, a study is done to the traversed path of an 
R-baryon in 1 m iron. The energies of R-hadrons which do not get stopped in the calorimeters is 
relevant, since those hadrons reach the outer parts of a typical detector; often the muon chambers. 
In Fig. 1111 the traversed length in iron is shown for different R-hadron masses, with a maximum 
length of 1 m. A rapid saturation takes place, i.e. above a rather moderate value R-hadrons will 
punch through a typical calorimeter. 
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Figure 9: Profile histograms representing the average R-hadron energy loss in one meter iron and 
the spread on this average as a function of its kinetic energy for the total energy loss (upper figures) 
and pure ionization losses for a singly charged R-hadron not undergoing nuclear interactions (lower 
figures). 

6 Signatures of heavy hadrons in ATLAS 

As an example of R-hadron signatures in a typical future detector, we consider the ATLAS detector. 
Concerning the selection of events containing R-hadrons, appropriate triggers would be the muon 
trigger, in case a charged R-hadron reaches them, or by the missing energy trigger, if not two R- 
hadrons are produced in a back-to-back topology in the transverse plane. Signatures would include 
the following. 

• Missing energy, the amount depending on the topology of the events. 

• A larger amount of ionization in the Transition Radiation Tracker, manifested by a large 
amount of high threshold hits, if the R-hadron moves slowly. This could mimic the transition 
radiation hits for particles with very large (3 values. 

• The E/p ratio, which is the amount of energy deposited in the calorimeters, divided by the 
momentum of the track, measured in the tracking system. This ratio is considerably smaller 
than that for light hadrons like pions, but larger than that for a muon, as is illustrated Fig. 1121 

• A characteristic longitudinal shower profile in the different components in the ATLAS calorime- 
ter, which would differ from that of pions, as well as from that of muons. 

• A characteristic transverse shower profile in the calorimeters. The low-energetic interactions 
result in a narrower shower than e.g. a shower caused by a high-energetic pion. 
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Figure 10: A comparison of total energy loss and the nuclear energy loss alone in 1 meter iron for 
an R-hadron with mass 100 GeV for a kinetic energy of 20 GeV (a) and of 300 GeV (b). The mean 
values of the distributions shown in Fig. (a) are 3.8 and 1.9 respectively, while those in Fig. (b) are 
13.8 and 12.9 respectively. 

• A large time-of- flight in the muon chambers. 

A detailed study of all the above aspects, as well as a detailed trigger efficiency study for the 
relevant triggers in the ATLAS experiment, will be discussed in a future publication. 



7 Conclusion 

In this paper, interactions of heavy hadrons in matter have been described, and their interactions 
have been simulated inside the GEANT 3 framework. Several approximations are made in the model 
with respect to cross sections, quark content, and interaction processes. The model presented here, 
as well as the simulation, should not be viewed as a final description, but provides a convenient 
platform on which to build in the future. 
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Figure 11: Comparison of traversed path of a singly charged R-baryon for different values of the 
R-hadron mass. The vertical lines represent the spread on the mean value. 
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Table 1: The mean value of the energy loss in 1 m iron for different values of Qq. The default value 
is Q = 1.1 GeV/c 2 



20 



P=100 GeV/c 

Muon 

R-hadron M=300 GeV/c" 

R-hadron M=20 GeV/c 2 

Pion 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 
E Id 

measured r measured 



Figure 12: The ratio E/p of the deposited energy in the ATLAS calorimeters and the momentum, 
measured in the ATLAS tracking system for particles generated with momentum p = 100 GeV/ 
The scale on the y-axis is arbitrary. 
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